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Abstract: We argue that the Weyl coordinates and the rod-structure employed to con- 
struct static axisymmetric solutions in higher dimensional Einstein gravity can be gen- 
eralized to the Einstein-Gauss-Bonnet theory. As a concrete application of the general 
formalism, we present numerical evidence for the existence of static black ring solutions 
in Einstein-Gauss-Bonnet theory in five spacetime dimensions. They approach asymptoti- 
cally the Minkowski background and are supported against collapse by a conical singularity 
in the form of a disk. An interesting feature of these solutions is that the Gauss-Bonnet 
term reduces the conical excess of the static black rings. Analogous to the Einstein-Gauss- 
Bonnet black strings, for a given mass the static black rings exist up to a maximal value 
of the Gauss-Bonnet coupling constant a'. Moreover, in the limit of large ring radius, the 
suitably rescaled black ring maximal value of a' and the black string maximal value of a' 
agree. 
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1. Introduction 



In recent years it has been realized that higher dimensions d > 4 allow for a rich landscape 
of black hole solutions that do not have four dimensional counterparts. The vacuum black 
ring solution of Emparan and Reall j^, [2| in d = 5 Einstein gravity is perhaps the best known 
example of such a configuration. The black ring has a horizon with topology S 2 x S 1 , while 
the Myers- Perry black hole |3| has a horizon topology S 3 . This solution provided also the 
first concrete piece of evidence that in higher dimensional gravity, the no-hair theorems 
of 3 + 1 dimensions do not apply. For example, in a 4 + 1 dimensional asymptotically 
flat spacetime with a given ADM mass and angular momentum, the geometry need not 
necessarily be that of the Myers-Perry black hole. 

The d = 5 Emparan and Reall black ring solution has been generalized in various 
directions, including configurations with abelian matter fields Physically interesting 

solutions describing superposed black objects (black saturns Q, bicycling black rings [10]. 
1 11], and concentric rings [jDl], [14]) were also constructed. However, in the static 
limit, all known d = 5 asymptotically flat solutions with a nonspherical horizon topology 
possess a conical singularity or other pathologies Q, Q, [||. 

All these results concern the case of d = 5 Einstein gravity theory and its various 
extensions with abelian matter fields. However, in five dimensions, the most general theory 
of gravity leading to second order field equations for the metric is the so-called Einstein- 
Gauss-Bonnet (EGB) theory, which contains quadratic powers of the curvature. The Gauss- 
Bonnet (GB) term appears as the first curvature stringy correction to general relativity [15 



16], when assuming that the tension of a string is large as compared to the energy scale of 



other variables. Inclusion of this term in the action leads to a variety of new features (see 
1 17], [18] for a review of the higher order gravity theories). 

Although the generalization of the spherically symmetric Schwarzschild solution in 
EGB theory has been known for quite a long time fl~9[] , the issue of axially symmetric 
solutions with a GB term is basically unexplored. In particular we do not know if the 
d > 5 black holes with a nonspherical topology of the horizon continue to exist when 
including stringy correction to the action. The main obstacle here seems to be that the 
Weyl formalism, which has proven so useful in the case of Einstein gravity, allowing for the 
discovery of a plethora of interesting exact solutions, has no straightforward extension in 
the presence of a GB term. 

In the absence of exact solutions, a natural way to approach this issue is to construct 
such configurations numerically. This paper aims at a first step in this direction, since we 
propose a framework for a special class of d = 5 static configurations with three commuting 
Killing vectors. In the absence of a GB term, this framework reduces to that used in Q to 
construct generalized Weyl solutions. Here we argue that some basic properties there are 
still valid in the presence of a GB term, in particular the rod structure of the solutions. 

As the simplest example of a d = 5 black object with a nonstandard topology of the 
event horizon, we present numerical evidence for the existence of static black rings in EGB 
theory. These solutions are found within a nonperturbative approach, by directly solving 
the second order field equations with suitable boundary conditions. These black rings 
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share most of the features of the Einstein gravity solution in Q. Although the inclusion 
of the GB term in the action reduces the conical excess, these configurations still possess 
an angular deficit. Moreover, for a given value of the mass, similar to the case of a d = 5 
EGB black string [20|, the static black ring solutions exist up to a maximal value of the 



GB coupling constant. 

The plan of the paper is the following. In the next Section we present a brief review 
of the Weyl formalism in Einstein gravity and argue that the coordinate system and the 
rod structure used there can be employed to construct EGB solutions as well. Section 
3 consists of a discussion of the Schwarzschild black hole and the uniform black string 
in EGB theory. There we present evidence that these configurations can also be viewed 
as generalized Weyl solutions within the framework of Section 2. Section 4 contains the 
main results of this work consisting of a systematic study of the static black rings in EGB 
theory. We give our conclusions and remarks in the final section. There we report also our 
results on charged generalizations of the static black rings present in EGB-Maxwell theory. 
In Appendix A we present some details on the EGB equations. Appendix B contains a 
discussion of some technical aspects involved in the numerical construction of the EGB 
static black rings. This includes a new coordinate system which has proven more suitable 
for the numerical study of the black ring solutions. 

2. The general formalism 

2.1 The Einstein-Gauss-Bonnet theory 

We consider the EGB action in five space-time dimensions 

where G is the five dimensional Newton constant and a' is the GB coefficient with dimen- 
sion (length) 2 . In string theory, the GB coefficient is positive, and this is the only case 
considered in this work 1 . R denotes the Ricci scalar and 

Lgb = R 2 - AR^RT + R^ pa R^ pa (2.2) 

the Gauss-Bonnet term with Ricci tensor R^ u and Riemann tensor R^pa. 

The variation of the action ( p.l| ) with respect to the metric tensor yields the EGB 
equations 

Ep U = Gp V + a'H^ = , (2.3) 



where 



G fit/ R^iv c^p^vR 



1 Also, a negative value of a' leads to a number of pathological features of the theory, see e.g. |l7[ . 
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For a well-defined variational principle, one has to supplement the action (^4]) with the 
Gibbons-Hawking surface term pl[ 

A E) = I d'x^K , (2.4) 



8vrG 



dM 



and its counterpart for Gauss-Bonnet gravity [16] 



7 > Cfl, =-4^.L d4 ^ j - 2G - K *'» • 



where ^y ab is the induced metric on the boundary, K is the trace of the extrinsic curvature 
of the boundary, G ab is the Einstein tensor of the metric j ab and J is the trace of the tensor 

J ab = 1 (2KK ac K c h + K cd K cd K ab - 2K ac K cd K db - K 2 K ab ) . (2.6) 

2.2 The Weyl solutions in d = 5 Einstein gravity and the rod structure 

Following the approach in @, we consider asymptotically flat, five-dimensional static and 
axisymmetric vacuum spacetimes with three commuting Killing vector fields Vu\ (i = 
1, 2, 3). The commutativity of Killing vectors [Vm, V(j\\ = enables us to find a coordinate 
system such that Vu\ = d/dx l and the metric is independent of the coordinates x 1 . In what 
follows, we shall put x 1 = t, x 2 = ip, and x 3 = (p. Then (d/dx 1 ) is the Killing vector field 
associated with time translation and (d/dx 2 ), (d/dx 5 ) denote the spacelike Killing vector 
fields with closed orbits. 

Here we invoke the particularization for d = 5 of the general theorem 2.1 in Ref. @: 
Let Vu-\,i = 1,2,3, be three- commuting Killing vector fields such that 

1. V [ { ^V^ ) Vg ) D u Vlj ) = holds at least at one point of the spacetime for a given i = 
1,2,3. 

2. V^R^V^VgVf^ 1 = holds for all i = 1, 2, 3. 

Then the two-planes orthogonal to the Killing vector fields Vu\,i = 1,2,3, are integrable. 

The first condition holds because we have assumed axisymmetry, while the second 
one is automatically satisfied as long as we restrict ourselves to the vacuum solutions of 
Einstein equations. 

As a result, the metric can be written in the canonical form [22] as 

ds 2 = e 2 ^ z \dp 2 + dz 2 ) + e 2U ^d^ 2 + e 2U3 ^dcp 2 - e 2U ^dt 2 , (2.7) 

where 0</?<oo, — oo < z < oo. Here it is most convenient to choose the three functions 
Ui as to satisfy the condition 

X>* = log p. (2.8) 

i 

This is compatible with the vacuum Einstein equations = (i = 1, 2, 3), which for the 
choice (^8) reduce to 

d 2 Ui 1 dUi d 2 Ui _ () 
dp 2 p dp dz 2 
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(the Einstein tensor for the metric ansatz (|2.7| ) is presented in Appendix A). One can 
see that ( |2.9| ) is just Laplace's equation in a (fictitious) three-dimensional flat space with 
metric ds 2 = dp 2 + p 2 d6 2 + dz 2 . 

From the other components of the Einstein equations Gp — G z z = and G z p = 0, 
we obtain the equations which determine the function v (p, z) for a given solution of the 
equation (|2.9|) 



v 



1 , P 



Y p + ^{U[ 2 + Ui 2 + U> 2 -U 2 -U> 2 -U 2 ), v = p{U' x + U' 2 + U' z ), (2.10) 



where a prime denotes the derivative with respect to p and a dot denotes the derivative 
with respect to z. Solutions with the ansatz Q2.7| ) and with U\,U2, U3 and v satisfying the 
equations (|2.9|), Q2.10 ) are usually called generalized Weyl solutions. 

Although the Einstein equations take a simple form in terms of (£7j, v), for the purposes 
of this paper it is more convenient to work with a set a functions /j defined as follows 

e 2 ^ = h(p,z), e 2U ^ =f 2 (p,z), e 2U ^ = f 3 (p,z), e 2U ^ = f (p, z). (2.11) 

This leads to a line element 

ds 2 = -f (p, z)dt 2 + h( P , z)(d P 2 + dz 2 ) + f 2 ( P , z)d^ 2 + f 3 (p, z)d<p 2 , (2.12) 

which was used in our study of the EGB static black ring solutions. 

In this paper we are mainly interested in configurations approaching asymptotically 
the five dimensional Minkowski spacetime, this being also the simplest solution of the 
equations (pOD , |2~lc| ). In this case, the metric functions fi have the following expression: 



fo(p,z) = l, Hp,z) = ^==, Hp,z) = Vp 2 + z 2 + z, Hp,z) = Vp 2 + z 2 -z. (2.13) 

2yV + z 2 

The usual form of the flat spacetime metric in the Hopf coordinates 

ds 2 = -dt 2 + dr 2 + r 2 {d6 2 + cos 2 6dtp 2 + sin 2 6dtp 2 ), (2.14) 
is found from ( |2.12| ), ( |2.13 ) via the coordinate transformation 

p= ^r 2 sin 20, z = ^r 2 cos29, (2.15) 

with 0<r<oo, 0<(9< vr/2. 

The equations (EO|), ( |2.10| ) possess a variety of physically interesting solutions. They 
can be uniquely characterized by the boundary conditions on the z— axis, known as the 
rod-structure pi, [22|, [23]. In pure Einstein gravity, the physically relevant solutions for Ui 
can also be thought of as Newtonian potentials produced by thin rods of zero thickness with 
linear mass density 1/2, placed on the axis of symmetry in the auxiliary three-dimensional 
flat space. Then the constraint fl2.8|) states that these sources must add up to give an 
infinite rod. 
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In this approach, the z— axis is divided into N intervals (called rods of the solution), 
[— oo, zi], [z\, zz],. ■ ■ , [zn-i, 00]. As proven in p2fl , in order to avoid curvature naked singu- 
larities at p = 0, it is a necessary condition that only one of the functions /o(0, z), /2(0, z), 
/3(0,z) becomes zero for a given rod, except for isolated points between the intervals. 

For the static case discussed here, a horizon corresponds to a timelike rod where 
/o(0, z) = while lim p ^o /o(p 5 z )/p 2 > 0. There are also spacelike rods corresponding to 
compact directions specified by the conditions f a (0,z) = 0, lim p _>o fa(p, z )/p 2 > 0, with 
a = 2, 3. A semi-infinite spacelike rod corresponds to an axis of rotation, the associated 
coordinate being a rotation angle. For example, the Minkowski spacetime ( 2.14| ) corre- 



sponds to two semi-infinite rods [—00, 0] and [0, 00]. Demanding regularity of the solutions 
at p = imposes a periodicity 2ir for both tjj and ip. (However, when several ip- or 93-rods 
are present, it may be impossible to satisfy simultaneously all the periodicity conditions). 

One of the main advantages of this approach is that the topology of the horizon is 
automatically imposed by the rod structure. This provides a simple way to construct a 
variety of solutions with nontrivial topology of the horizon (including multi-black objects). 
Since ( |2.9| ) is linear, one can superpose different solutions for the same potential The 
nonlinear nature of the Einstein gravity manifests itself through the equation ( 2.1 0| ) for the 
metric function v. 

2.3 The static, axially symmetric Einstein-Gauss-Bonnet configurations 
2.3.1 The equations 

One of the main purposes of this work is argue that the ansatz ( |2.7| ) and the associated rod 
structure can be used to construct physically relevant solutions in EGB theory. A priori, 
it is not clear that this metric ansatz is valid also in this case, since the second assumption 
in the general theorem mentioned above {i.e. ^V^V^Y^ = 0) does not hold in 
general for EGB theory. Thus the situation here is similar to the case of caged black holes 



1 24], or nonuniform black strings [26|, [27] in Kaluza-Klein Einstein gravity, where 



the validity of the metric ansatz could be proven only a posteriori, after solving the field 
equations. 

The equations for the functions f\, f^-, an d /o are found by using a suitable combi- 
nation of the EGB equations, E\ = 0, E p p + E z z = 0, E't = 0, and E$ = 0. Details on these 
equations and the explicit form of the tensors and are presented in Appendix A. 

The remaining equations E§ = 0, E p p — E* = yield two constraints. Following Pq] , 
we note that setting E\ = E* = E p p + E z z = in V p E^ = and V p E» z = 0, we obtain 
the Cauchy-Riemann relations 

d p (V=P£) + d p (^g\{EP - e%)\ = 0, d p (V=P£) - d z (V=4(££ - K)) = 0. (2.16) 

Thus the weighted constraints satisfy Laplace equations, and the constraints are fulfilled, 



when one of them is satisfied on the boundary and the other at a single point [26|. 

Due to the GB contributions, the second order equations for the functions /, are much 
more complicated than in the case of Einstein gravity and do not reduce to the simple 
Laplace equation. Thus the functions Ui are no longer harmonic. Also, one can verify 
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that the central property (|2.8| ) of the Einstein gravity Weyl-solutions does not hold in the 
presence of a GB term and thus one cannot set /0/2/3 = p 2 - 

Therefore finding closed form solutions within this approach looks unlikely. However, 
the solutions can be constructed numerically, by solving boundary value problems. A major 
advantage of the ansatz ( |2.7| ) is that the (p, z) coordinates have a rectangular boundary in 
which all boundaries coincide with the coordinate lines and thus are suitable for numerics 2 . 

2.3.2 The rod structure 

The central point in this approach is that the rod structure, as explained above for the 
case of Einstein gravity, can be used also for solutions of the EGB theory 3 . This would fix 
the boundary conditions along the z— axis for the functions and thus the topology of the 
horizon. 

Here one starts by noticing that the following generic form of a solution near the z— axis 
is compatible with the EGB equations: 

f l (p,z)=f i0 (z)+p 2 f i2 (z) + ..., (2.17) 

where the functions fik(z) are solutions of a complicated set of nonlinear second order 
ordinary differential equations. Then, similar to the case of Einstein gravity, the z— axis is 
divided into ./V intervals-the rods of the solution. Except for isolated points between the 
rods, one assumes that only one of the functions /o(0, z), f 2 (0, z), f 3 (0,z) becomes zero 
for a given rod, while the remaining functions stay finite at p = in general. (In fact, if 
more than one of these functions is going to zero for a given z inside a rod, one can prove 



following the arguments in [22], that there is a curvature singularity at that point.) Again, 
one imposes the condition that the ./V intervals must add up to give an infinite rod. 

For example, for a rod in the ^-direction, one finds the following expansion of the 
metric functions as p —* 0: 

Up, z) = f 20 (z) + p 2 f 22 (z) + ..., A(p,z) = f 10 (z) + p 2 f 12 (z) + ..., (2.18) 
f 2 {p, z) = p 2 f 22 (z) + p 4 f M (z) + ..., f 3 (p, z) = f 30 (z) + p 2 f 32 (z) + . . . . 

The important feature here is that the constraint equation E z p = implies fio(z)/f 22 (z) = 
ci, i.e. a well-defined periodicity for the coordinate ip. 

Therefore, in order to cure the conical singularity at the rod, the coordinate ip should 
have a periodicity Aip = 2ir^/~c~[. A different periodicity of tp implies the occurrence of a 
conical singularity (this is the case if there are are several different ip-io&s). 

Similar results holds also for a rod in the (^-direction (i.e. when interchanging f 2 and 
^3), the periodicity of 92 there being again fixed by the constraint equation E z p = 0, i.e. 
lim p ^ /O 2 /i//3 = c 2 . 



2 This is not the case of the ring coordinates used in most of the studies on black ring solutions. The 
spatial infinity corresponds there to a single point. 

^However, note that the interpretation of a rod as corresponding to a zero thickness source with linear 
mass density 1/2, placed on the axis of symmetry in a auxiliary three-dimensional flat space is no longer 
valid in EGB theory. 
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Schwarzschild- Gauss- Bonnet black hole black ring 

(a) (b) 



black string black Saturn 

(c) (d) 



Figure 1. The rod structure of the solutions is shown for several EGB solutions. The thin lines 
denote the z— axis and the thick lines denote the rods. 

A finite timelike rod corresponds to an event horizon, where 4 

f (p, z) = p 2 f 02 (z) + p 4 f 04 (z) + ..., A(p,z) = f 10 (z) + p 2 f l2 (z) + ..., (2.19) 
/ 2 (p, z) = f ^(z) + p 2 f 22 (z) + ..., h(p,z)= f 30 (z) + p 2 h 2 {z) + ..., 

with linip^o P 2 fi/fo = c 3i which fixes the Hawking temperature of solutions. 

Thus, depending on the physical situation we consider, the boundary conditions along 
the z— axis are fixed by the above relations. The obvious boundary conditions for large p, z 
is that the functions fa approach the Minkowski background functions ( 2,13| ). 



Similar to the case of Einstein gravity, the topology of the horizon is fixed by the 
boundary conditions satisfied by f 2 and at the ends of the corresponding (finite) timelike 
rod 5 . 

For example, if either end of this rod continues with rods of different angular directions, 
then the event horizon has an S 13 topology (see Figure la). A black ring corresponds to 
fi ° r vanishing at both ends of the finite timelike rod associated with the horizon (see 
Figure lb). One can consider as well a black Saturn combining both types of black objects 



4 fik(z) here should not be confused with those in (2.18). 

5 A timelike rod extending to infinity corresponds to an acceleration horizon. 
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above, with two different horizons (see Figure Id). Moreover, if both ip— and ip— rods 
extend to infinity, then the spacetime is asymptotically flat. Solutions in Kaluza-Klein 
theory (for example a black string without a rod on the ip— or direction (see Figure 
lc)) can also be considered. 

It is tempting to conjecture that, similar to the case of Einstein gravity [23], a d = 5 
EGB solution within the ansatz ( |2.12| ), is uniquely specified by its rod structure. 

2.3.3 The physical quantities 



The general results in the literature [28] show that, similar to the case of Einstein gravity, 
the mass M of an asymptotically flat EGB solution can be read from the asymptotic 
expression for the metric component gu 

\GM 

-<to = /o~l-— -== + ... ■ (2.20) 
zn\J + z l 

Supposing we have an event horizon for z\ < z < Z2, the horizon metric is given by 6 

da 2 = A(0, z)dz 2 + / 2 (0, z)dip 2 + / 3 (0, z) dip 2 . (2.21) 

Two quantities associated with the event horizon are the event horizon area Ah and the 
Hawking temperature. For the metric ansatz ( 2.12 ) these are given by 



A H = AipAip r dz v / f 1 (0,z)f 2 (0,z)f 3 (0,z), Th = -— lim J {°} P / Z \ . (2.22) 
For solutions in EGB theory it is also of interest to compute the Ricci scalar of the horizon 

o . . 1 ( /i(0,z)/2(0,z) , fi(0,z)f 3 (0,z) _ f 2 (0,z)f 3 (0,z) fi(0,z) /|(0,z) _ o / / 2 (0,z) M0,z) s \ (0 9 oA 



2/l(0,z) yfi{0,z)h(P,z) ~ r h(0,z)f 3 (0,z) f 2 {0,z)f 3 (0,z) ^ J^z) ^ JJ^z) ^ f 2 (0,z) ~ r h(0,z) 

Considering now the case of a space-like ip— rod for z% < z < 24, one writes the line element 
on this three-dimensional surface £5 

da 2 = /i(0, z)dz 2 + / 3 (0, z)dp 2 - / (0, z)dt 2 . (2.24) 

The first quantity of interest is the proper length of the rod 

rz4 

L= / dzVhKzj, (2.25) 

J z 3 

(note that for a finite rod, L differs from the coordinate distance Az = z^ — 23). 

The solutions we are interested in may possess a conical singularity along some region 
of the symmetry axis. To define a conical singularity for a rotational axis with angle ip one 
computes the proper circumference C around the axis and its proper radius R and defines: 



dC 

OL = — ^7 

dR 



= lim = lim ^v^^ ( 2 . 26 ) 

R =o Jo y/gT P dp p^o ^/gyp 



°If there are several horizons, then one should write such an expansion for each of them. 
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where Aip is the period of tp. The presence of a conical singularity is now expressed 7 by 
means of: 



S = 2tt - a = 2tt f 1 - lim J f 2 2 j P / Z \ ) , (2.27) 

such that 5 > corresponds to a conical deficit, while 5 < corresponds to a conical excess. 
A conical deficit can be interpreted as a string stretched along on a certain segment of the 
z— axis, while a conical excess is a strut pushing apart the rods connected to that segment 
(in fact, for d = 5, the struts and strings are two dimensional surfaces). Similar to Einstein 
gravity, a constant rescalings of ip can be used to eliminate possible conical singularities on 
a given segment, but in general, once this is fixed, there will remain conical singularities 
at other ^-segments. 

For 5 < 0, we have found it convenient to introduce the quantity 

which has a finite range and measures the 'relative angular excess'. 

Of interest here is to compute the spacetime area spanned by the ^-rod. This is done 
by going to the Euclidean section t —* it and evaluating the quantity 

Area = pAip / dsv7o(0,*)/i(0, z)f 3 (0, z), (2.29) 

where (3 = 1/Tfj is the periodicity of the Euclidean time. For completeness, we give here 
also the expression for the Ricci scalar on the ip— rod, 

r> . 1 ( MMMM _l A (0,^/3(0^) _ /o(0,z)/3(0,*) , ft(0,z) , fi(0,z) _ q/ /o(0,*) , iaCM s \ , 

— 2fi!Q.z) fi(Q.z)f n (Q.z) fUQ.z)h(Q.z) fn(Q.z) MO.z) ^ TMTTl ^ f?(0 ^ Z[ > fn(Q.z) ^ MQ.z) > ' V z - OL V 



S 2/i (0,«) yh(0,z)f (0,z) ^ fi(0,z)f 3 (0,z) /o(0,z)/ 3 (0,«) ^ ffi&z) ^ 7f(M ^/<>(0,z) ^ f 3 (0,z) 

Of course, similar expressions can be written when considering instead a (p-rod. 
2.3.4 Remarks on the free energy and thermodynamics 

The discussion in this subsection applies to asymptotically flat solutions (although it can 
easily be generalized to the Kaluza-Klein case). 

The gravitational thermodynamics of the EGB black objects can be formulated via 
the path integral approach [21, 29|. In what follows it is important to use the observation 
that one can write 




/2/3 



R\^-g = -- I W^/o)+ W^jf /o) 1 . ( 2 - 31 ) 



{Hi + \L GB )^-g = - {8 p T p + ihT, i . (2.32) 



7 Note that, in some sence, fixing 8 is the analogue of computing the Hawking temperature on the 
Euclidean section. 
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where 



rp _ j foh ( /0/ /2/3 hh ^ , fo f hfo , /3/0 , ^2/3 , Tg/g N , fo_ tfif/j_ff^ 
7o \/i 72/3 72/3 7i 7o72 7o73 7273 72/3 7i72 

+ 757" (/2/1 ~ /1/2) + -prrU3.fi ~ /1/3) + 72V (/1/3 + A/3) 
7i 72 7i 73 7i 73 

7^(| + |) " + + iM + |) ) , C2.:J:il 

7o7i 73 72 7i 72 73 7i 72 73 / 

hh ( fo / i f i , i f i\ , /o /; ; ww\ , /0/0//3 , /2> 



^ = V — T~ ( TTT^f'z + 73/2) + 7-7-^(72/3 - 72/3) + 7"7~(7" + 7") 
V /o V/1/2/3 /1/2/3 /0/1 73 72 

~~ T7~(7^ + T 1 ) + 727"(7i/2 _ /2/1) + 72V (/1/2 + A/2) 
7o/i /2 73 7l/2 7l/2 

+ t4(A£ - A70 + + A/3) - ^(f + 1) + ^4 + r)| 

7i/3 7i/3 7i 72 73 7o 72 73 / 

When computing the classical bulk action evaluated on the equations of motion, one re- 
places the 72 + a 1 Lqb volume term with 2(72* + a'(Hj + Lqb/2)) and make use of ( 2.31| ) 



to express it as a difference of two boundary integrals. The boundary integral at infin- 
ity should be evaluated together with the contributions from 7^ and 7^ .As usual, 
this quantity is divergent. To regularize it, one has to subtract the contribution of the 
Minkowski background for both Einstein and Gauss-Bonnet boundary terms 8 . 

A direct computation implies the following expression for the tree level action of a single 
black object in d = 5 EGB theory (the extension to multi-black objects is straightforward) 

7 = 0(M - ^T H (A H + c!A x )), (2.34) 

with 

j3 



Ax = 2 drxVhRx h , (2.35) 

where h = V/1/2/3 is the determinant of the induced metric on the horizon and R^ h is 
the event horizon curvature as given by ( |2.23| ). 

The above results hold for the case of configurations with a regular z— axis (i.e. the 
periodicity of both angles ip and ip is 2ir everywhere). It is interesting to extend this 
analysis to the case of EGB solutions with conical singularities (this is the case of the 
black rings discussed in the Section 4). For simplicity, we shall consider here the case of a 



Note that for asymptotically flat solutions, the Gibbons-Hawking boundary term (2.4) gives a divergent 
contribution to Jo, while the GB boundary term reduces to a constant factor. It would be interesting to 
generalize the quasilocal formalism and the renormalized boundary stress-tensor from Einstein gravity to 
EGB theory. This will avoid the requirement to choose a background for the solutions. For example, we 
have found that the usual counterterm used in |$(J] for the black rings in Einstein theory regularizes also 
the mass and action of some (asymptotically flat-) EGB solutions, including the black rings discussed in 
Section 4. However, this approach implies the existence of a constant term ~ /3a' in the action, originating 
in the contribution of the GB boundary term (Ejj). 
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single singular section of the z— axis associated with a ip— rod (the generalization to other 
cases is straightforward). Then the conical singularity will add an extra contribution to 
the total tree level Euclidean action of the system, which leads to a more complicated 
thermodynamics of the system (see e.g. [31] for a related discussion for the d = 4 Israel- 



Kahn solution). This contribution can be evaluated by using the relations [32] 

- / d 5 Xy/gR = Area 5, - [ d 5 x^/gL GB = Area! 5, (2.36) 
where Area is the space-time area of the surface spanned by the conical singularity and 

Area\ = 2 / d?x VhRxs, (2.37) 



with h = V /1/0/3 the determinant of the induced metric on the surface T,$ and Ry, s the 
Ricci scalar on as given by ( 2.3C| ). 



Thus the total action in the presence of a conical singularity becomes 

I = 0(M - -j—Tji(Ah + a'Ai)) - -^—{Area + a' ' Area x ) 5. (2.38) 

The free energy of the solutions is identified as F = ThI = M. — ThS, where M. is the mass 
which enters the thermodynamics and S is the entropy. The entropy of the EGB black hole 
solutions without conical singularities can be written as an integral over the event horizon 



S = -^f d 3 x^/I(l + 2a'Rx h ), (2.39) 

which is the sum of one quarter of the event horizon area plus a Gauss-Bonnet correction. 
In this case, the mass M computed from the first law of thermodynamics is equal to the 
mass M computed at infinity. 

It would be interesting to perform a similar computation for solutions with conical 
singularities. Similar to the case of Einstein gravity |34|], a new extensive parameter A = 
{Area + a' Area\) / (3 associated with the rod containing the conical singularity will appears 
here, while the free energy becomes a function of both Tjj and A. Also, the first law 
of thermodynamics contains an extra- work term TdA (with T = —5/8ttG the tension 
associated with A) ■ Then the entropy S and the thermodynamical mass A4 of the physical 
system are given by 



s dF 



dT H 



, M = F + T H S. (2.40) 

A 



The thermodynamical stability of the EGB solutions can be studied in the usual way. For 
example, one defines the specific heat of a black object C = Th ^^^"^j with C > for 
thermodynamically stable solutions. 



- 12 - 



3. The known static solutions in d = 5 EGB theory 



3.1 The Schwarzschild black hole in EGB theory 
3.1.1 The solution in Schwarzschild coordinates 

The black hole solutions of EGB gravity have been studied begining with the work of 
more than 20 years ago. Most of the solutions in the literature restricted to static, 
spherically symmetric configurations and a Schwarzschild coordinate system. A suitable 
metric ansatz in this case is 

dr 2 

ds 2 = — - + r 2 (d6 2 + sin 2 6dip 2 + cos 2 0# 2 ) - N(r)a 2 (r)dt 2 , (3.1) 
N(r) 

and the expressions for N(r) and a(r) depend on the matter content of the theory. 

The Hawking temperature of a generic black hole ( [3.1] ) is Tjj = A^'(r/ l )a(r/ l )/(47r), 
where a prime denotes the derivative with respect to the radial coordinate and is the 
largest positive root of N(r), typically associated to the outer horizon of a black hole. 
Also, the event horizon area of a black hole is Ah = V^r^ (with V3 = 2-7T 2 the area of the 
three-sphere) . 

The complexity of the EGB theory, basically due to higher order terms in the curvature 
tensor, makes the task of finding exact solutions very difficult. To the best of our knowledge, 
the only d = 5 static, spherically symmetric EGB solutions known in closed form are the 
generalizations of (electro-) vacuum Einstein gravity configurations 9 . The Schwarzschild 
black hole in EGB theory has o(r) = 1, the metric function N(r) being given by 10 



r 



2 / a^//^'2 



8a' {r' 2 + 2a') 



Ar M" 1 + M 1 ~V 1+ f I ' (3 ' 2) 

The constant in ( |3.2| ) corresponds to the event horizon radius. One should note that 
the EGB Schwarzschild solution exists for all > and a' > 0. As r —* one finds 



+ 4a' 



At short distances, the EGB Schwarzschild solution and its Einstein gravity counter- 
part are substantially different due to the effect of the GB term. An interesting fea- 
ture is that the EGB Schwarzschild metric turns out to be finite at the origin, since 
N(r) — ► 1 — J (r 2 , + 2a') /(2a') as r — > 0. However, at large distances (r 2 S> a') the 
EGB black hole behaves like the Schwarzschild solution. Once the event horizon radius 
is fixed (as is done in the numerical construction of the solutions), the parameter a' enters 
the 1/r 2 term in the asymptotics, 



9 Exact solutions describing cosmic strings in d — 5 EGB theory were found in pfj]. Obviously, these 
solutions cannot be described within the metric ansatz (3.1). 



The usual form for N(r) in the literature is in terms of r h = m — 2a', which corresponds to fixing the 



mass of the solutions. However, the expression (3.2) fits better with the purposes of this work 
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Also, as a' — > 0, one recovers to leading order the Schwarzschild black hole expression 

r 2 2 2r 4 

N(r) = 1 - -A + (-- + -Ay + O(a') 2 - (3.5) 
The Hawking temperature and the mass of these solutions are given by 

The above expression for M shows the existence of a mass gap: black holes exist only for 
M > SV^a' /(8nG). Thus the mass spectrum of the spherically symmetric EGB black holes 
is bounded from below. 

A straightforward computation leads to the following expression for the entropy of the 
EGB Schwarzschild black holes 

S = S + S C with S = ^, S c = 3a'^. (3.7) 

One can easily verify that the first law of thermodynamics dM = TndS also holds. 

Without entering into details, we mention the existence of some substantial differences 
between the thermodynamics of the EGB Schwarzschild solutions and their Einstein gravity 
counterparts. If the black holes are large enough 3> \fa! ', then they behave like their 
Schwarzschild- Tangerlini counterparts. A different picture is found for small values of r^, 
since Tjj — rh/(8ira') in that case. Therefore, the specific heat changes its sign at length 
scales of order ~ \[a! . This implies the existence of a branch of five dimensional EGB 
black holes which is thermodynamically stable (see the Ref. [17| for a review of these 
aspects). 

3.1.2 The solution in Weyl-type coordinates 



To bring the generic metric fl3.1|) to Weyl form (2.12), one considers the coordinate trans- 
formation 

p = -rjj sinh G(r) sin 26, z = -rjj cosh G(r) cos 26, (3-8) 

where ro is defined by the asymptotic expansion of N(r) = 1 — (ro/r) 2 + . . . (i.e. Tq = 
r\ + 2a'). The function G(r) above is defined as 

G(r) = 2 ^ _*L= (3.9) 

Jr h XyJN(x) 

(with G(rh) = 0). Then the Schwarzschild coordinate r is expressed in terms of Weyl 
coordinates p, z as 

r(p,z) = G~ l (arsinh^X(p, z)), (3.10) 

where 



xm = i V + + JX + (i - l -^^f) • (3.n; 



2 
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Figure 2. The metric functions of an EGB Schwarzschild black hole are shown versus z for 
several values of p. The relevant parameters here are a = 1, a' = 0.01. 



A straightforward but cumbersome computation leads to the following expression for the 
metric functions for the parametrization (2.12) 



f (p,z) = F(r(p,z))a 2 (r(p,z)), fi(p,z) 



r(p,z)P 2 (p,z) 
1 + P*(p,z) ' 



(3.12) 



f 2 (p,z) = ^r 2 (p,z)(l + 1 = 



), hip, z) = \r 2 {p, z)(l 1 == ), 



with 



P(j>,z) 



p X(p,z) + 1 
z\ X(p,z) ' 



(3.13) 



Unfortunately, the integral (|3.9| ) cannot be computed in closed form for the known solutions 
of EGB theory (except for the asymptotic expressions as r — * and r — * oo). However, 
the expressions (|3.9|) , ( |3.12j ) can be evaluated numerically. 
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For any value of a', the rod structure of the EGB Schwarzschild black hole as resulting 
from ( 3.12| ) consists of a semi-infinite space-like rod [—00,— a] (with /2(C), z) = there), 



a finite time-like rod [—a, a] (/o(0, z) = 0) and a semi-infinite space-like rod [a, 00] (with 
vanishing /3(C), z)) in the ^-direction (with a = Tq/4). Thus the topology of the horizon 
is S 3 as required (see Figure la). A plot of the metric functions /, exhibiting this rod 
structure for a typical EGB Schwarzschild solution is shown in Figure 2. In principle, 
most of the physically relevant properties of the EGB Schwarzschild black hole can also 
be rederived within the metric ansatz ( |2.12| ), However, the required computation is much 
more difficult for that coordinate system. 

3.2 The static uniform black string 

It is of interest to briefly review the situation for a different type of black object in EGB 
theory which can also be studied within the ansatz (2.12). In Einstein gravity, one can 



construct uniform d = 5 black string solutions by adding a flat direction to any d = 4 
vacuum black hole. (These solutions can still be written in the d = 5 Weyl form ( 2,12; ), 



the new direction tp being trivial.) However, it is straightforward to check that this simple 
construction does not work in the presence of a GB term in the action 11 , unless the solutions 
are conformally flat. Although no exact solutions describing d = 5 black strings in EGB 



theory are known so far, these configurations were studied numerically in Ref. [2C]. These 
solutions can be constructed within a metric ansatz related to fl3.1| ) 

ds 2 = — - + r 2 (d9 2 + sin 2 6d^ 2 ) - N(r)a 2 (r)dt 2 + b(r)dip 2 , (3.14) 
N(r) 

(note that < 9 < tt in this case while the periodicity of tp is not fixed a priori). The new 
feature here as compared to the case of Einstein gravity is that the metric component Qipip 
differs from one. 

The event horizon of a uniform black string is located at r = r/, > 0, where the 
following approximate form of the metric functions holds 

N(r) =N 1 {r-r h )+0(r-r h ) 2 , a{r) = a h + a x (r - r h ) + 0(r - r h ) 2 , (3.15) 
b(r) =b h + 61 (r - r h ) + 0(r - 

where 

_ 2b h (r h - N ir 2 + 4N ia ' + 2N 2 r h a' 



r h ) 2 



N 1 (r 2 +Aa')(2N 1 a' - r h ) 
ax =aJ- r 2 h {-l + N irh ){l + Wir h ) + 2N x r h {A + 3^(3 + 4N irh ))a' (3.16) 

- 8^(1 + 3A^ir h (-3 + N irh ))a' 2 - lUNfa' 3 ^ (w?r h (rl + Aa'){-2a' + {r h - 2N ia ') 2 )^ 



-1 



while 



r h ~ \ r h~ 8a ' 
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Figure 3. The metric functions fi of a typical uniform black string solution in EGB theory 
are shown versus z for several values of p. The relevant parameters here are a = 1, a' = 0.175. 

The free parameters in the near horizon expansion are bh > and > 0. 

The above relations imply the existence of a minimum horizon size for a given value 
of the GB coupling constant 



(3.18) 



Since, as discussed in Ref. pCp , the horizon radius is decreasing monotonically with the 
mass of the solutions, the relation ( |3.18| ) shows again the existence of a minimal value of 
the mass for a given GB coupling constant a' . 

At infinity, the background approached by a black string is the four-dimensional space- 
time times the cp— direction, ds 2 = dr 2 + r 2 {d0 2 + sin 2 9dij) 2 ) — dt 2 + dip 2 . The solution as 



It is interesting to notice that this is valid also for solutions with a cosmological constant. 



r — ► oo is written in terms of two parameters Q, c 2 : 

N(r) = 1 + C -1^L + + (l/r% a(r) = 1 - °f + ^ " ° t] + 0(l/r 3 ), 
r 4r 2r or 

6(r) = l + ^ + ^i + 0(l/r 3 ). (3.19) 

Similar to the case of Einstein gravity, the EGB black string solutions possess two 
global charges - the mass M and the tension T, associated with the Killing vectors d/dt 
and d/dip, respectively. These global charges are fixed by the constants q, c z in the 
asymptotic expansion ( 3.19| ): 



M = ^ {2Ct ~ C ^ r= 4^G (Q - 2c ^ (3 - 20) 

Without entering into details, we mention that the black string solutions can also be 
recovered within the metric ansatz ( 2.12| ). The coordinate transformation between this 



ansatz and ( 3.14 ) can be worked out in a similar way to (|3.8| )-( 3~T^ ) (note however that 



p = ^ro sinhG(r) sin 9, and z = cosh G(r) cos 9 in this case). These solutions have a 
semi-infinite space-like rod [— oo,— a] along the vp— direction, a finite time-like rod [—a, a] 
corrresponding to the event horizon and a second semi-infinite space-like rod [a, oo] again 
in the ^-direction (thus there is no rod along the (p— direction, see Figure lc). This can 
also be seen in Figure 3, where the metric functions of a tyical EGB uniform black string 
solution are shown as a function of z for several values of p (one can notice the nontrivial 
shape of the metric function f^). 

However, the EGB black string solutions are much more difficult to study in the co- 
ordinate system ( p. 12 ), since in this case one deals with partial differential equations. In 



particular, it is much more difficult to prove analytically the existence of a maximal value 
of a' for a given length of the finite timelike rod 12 . 



4. The static EGB black rings 

The physical intuition (supported by the results in the previous Section) suggests that 
all known solutions in Einstein gravity admit generalizations in EGB theory. While it is 
rather inconvenient to use the metric ansatz Q2.12 ) for the study of EGB Schwarzschild 



black holes and EGB black strings, this is not the case for more complicated solutions 
with a nonspherical topology of the horizon. In fact, in our opinion, the simplest way 
to construct EGB generalizations of such objects is within the metric ansatz ( [2.12| ), by 
imposing the same rod structure as in the absence of the GB term. 

In this Section we present numerical evidence for the existence of static black rings, as 
the simplest example of an asymptotically flat black object with a nonspherical topology 
of the horizon in EGB theory. 



12 However, the coordinate system ( |2.12| ) makes possible to attempt a numerical construction of more 
complicated Kaluza-Klein solutions in EGB theory, e.g. multi-black strings or configurations with bubbles. 
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4.1 The static black rings in Einstein gravity 

The rod structure of a static black ring solution in Einstein gravity is exhibited in Figure lb. 
It consists of a semi-infinite space-like rod [—00,21] in the ^-direction (thus /2(C), z) = 
there), a finite time-like rod [^1,^2] (/o(O>*0 = 0), a second (and finite) space- like rod 
[2:2,2:3] in the ^-direction, where /2(0, z) = again, and a semi-infinite space-like rod 
[23,00] (/3(C), z) = 0) in the (^-direction (and z\ < Z2 < z 3 ). 

The metric functions /, of the static black ring 13 are given by f2|,[f2^] 

, #2 + 6 , (R l + £ 1 +R 2 -&)(0--c)R 1 + (l + c)R2 + 2cR 3 ) 

JO — ~Fi — 7~Ti /l 



Ri+ii 8(1 + c)R 1 R 2 R 3 

-6) (Rs- 
Ri-ti 



h= ^~t m+C3 \ A = (4.1) 



where 



Ci = z - Zi, Ri = \ p 2 + Q and z\ = -a, z 2 = a, z 3 = b, (4.2) 



a and b being two positive constants, with c = a/b < 1. Roughly speaking, a fixes the size 
of the horizon, while b provides a measure of the radius of the ring's S 1 . 

Since the orbits of tp shrink to zero at —a and a while those of tp do not vanish anywhere 
there, the topology of the horizon is 5 2 x 5 1 , see Figure lb (although the S 2 is distorted 
away from perfect sphericity). 

The mass, event horizon area and Hawking temperature of this solution are: 



Although the static black ring solution is asymptotically flat 14 , it contains a conical singu- 
larity for the finite t/'-rod, since 5 as defined by ( |2.27 ) is nonzero: 



One can easily see that this is a negative quantity 5 < 0, which implies the existence 
of a two-dimensional disk-like deficit membrane (with negative deficit) that prevent the 
configuration from collapsing. 

Another quantity of interest is the area of the spacetime spanned by the conical singu- 
larity and the proper length of the finite ip-rod, which are computed according to ( 2.2S| ), 
d2T25|): 



Ar P n = R 27T ( b Q ) 3 /2 L = ^2\fb~^a E(n), (4.5) 
ya + b 



13 Note that the function /i(0, z) behaves as l/\z — z%\ as z — > Zi. 

14 The Einstein gravity static black ring solution in Q admits an alternative interpretation as a ring 
sitting on the rim of a membrane that extends to infinity. (This is found by requiring that the periodicity 
of tp is 2tt on the finite tp-rod.) However, the asymptotic metric is a deficit membrane in this case. 
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where n = (b — a)/(b + a), E(n) being the complete elliptic integral of the second kind. 

In terms of the dimensionless parameter a/b, one may think of a static black ring 
as interpolating between two limits. As a/b — ► 1, the finite ip— rod vanishes and the 
Schwarzschild metric is approached, with 5 — > — oo. As a/b — > 0, the second ip— rod 
extends to infinity and the solution becomes, after a suitable rescaling 15 , a black string, 
i.e. the four dimensional Schwarzschild black hole uplifted to five dimensions. 

4.2 The static black rings in EGB theory 

4.2.1 The ansatz and the physically relevant quantities 

The EGB generalizations of the Emparan-Reall black rings are found by solving the EGB 



equations for the metric ansatz ( 2,12| ) 



The boundary conditions satisfied by the EGB black ring metric functions are similar 
to those in Einstein gravity. At p = 0, the function /o vanishes for —a < z < a (i.e. on 
the horizon), fi is zero for — oo < z < —a and a < z < b, while f% vanishes for z > b. As 
a result, along the horizon the orbits of tp shrink to zero at z = —a and z = a, while the 
orbits of cp do not shrink to zero anywhere. Thus the topology of the horizon is S 2 x S . 
From ( 2.1£| ), ( |2.19| ), one can see that, for a given rod with one of the functions vanishing, 



f a = 0, the other /, satisfy Neumann-type boundary conditions, d p fi\ p= o = (with i ^ a). 
At infinity, we require that the functions fi approach the Minkowski form ( p,13[ ). 
In practice, we have found it convenient to take 

fi = fiFi, (4-6) 

where ff are background functions, given by the metric functions of the Einstein gravity 
black ring solution ( |4.1[) . The advantage of this approach is that the fi will automatically 
satisfy the desired rod structure. Moreover, this choice 'absorbes' the divergencies of the 
functions /2 and fs as r — » oo coming from the imposed asymptotic behaviour. 

The equations satisfied by the Fi can easily be derived from the general set of EGB 
equations 16 . As for the boundary conditions, the relations ( [2.18D , ( |2.19| ) together with the 
expressions (|4.1| ) of the background functions ff imply 

d p Fi\ p= Q = 0, for — oo < z < oo, 

and Fi = 1 as p — > oo or z — ► ±oo. 

The constraint equation E z p = results in F<i/F\ = const, on the ip-rods. Now, 
to be consistent with the assumption of asymptotic flatness, one finds const. = 1 for 
— oo < z < —a. The value of this ratio for the second rod with a < z < b is obtained only 
as a result of the numerical solution. A similar reasoning implies F^/F\ = 1 on the 99-rod 
(b < z < 00). 



15 For the line element used in this work, this rescaling is r — > v2&f, z —* V2bz, <p — > tp/v2b together 
with a — > V2ba. 

16 Note, however, that the field equations become much more complicated in terms of the Fi, with the 
number of terms increasing drastically. 
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All of the physically relevant quantities except for the mass are encoded in the values 
of the functions Fi, ff at p = 0. The event horizon area of static black rings in EGB theory 
is given by 



A H = A^2n dzy/fthh = AtfjiTra^^ J* dz y/(b - z)F 1 F 2 F 3 , (4.7) 

where Aip is the periodicity of the angular coordinate ip on the horizon. 

The Hawking temperature can be computed from the surface gravity or by requiring 
regularity on the Euclidean section 



1 a + b F 

T " = 4^V^-ViV (48) 

where the constraint equation E z p guarantees that the ratio Fq/F\ is constant on the event 
horizon. 

At infinity, the five dimensional Minkowski background is approached, with Aip = 2ir 
there. From fl2.2CQ , the mass M of the solutions can be read from the subleading term in 
the asymptotic expression for /o = /qFq. To obtain a Smarr-like relation (which is useful 
in numerics, see Appendix B.2), we consider the EGB equation E\ = G\ + a' Hi = in the 
form R\ = 1/2R — a'Hj and integrate over a spacelike hypersurface. The volume integral 
on the left hand side of the equation reduces to surface integrals at the horizon and at 
infinity, which can be evaluated in terms of the mass, area and Hawking temperature, 

J RlV=gd 4 x = X - hnTnAa - yGAfj . (4.9) 

Substituting this expression on the left hand side and solving for the mass yields 

mOM = m + l a , = m -SJ{R- 2«^}» (4.10) 

In the limit a' — > the integral l a t vanishes and the relation reduces to the usual Smarr 
relation. 

To further characterize the properties of the horizon, we introduce the minimal and 
maximal S 1 horizon radii, i? m i n and -Rmax; defined via 

-27T 



Rutin = y j V9^\ P =o,z=a d V = V^b-a)F 3 (0,a), (4.11) 

J 

l - j V9^\ P=0 , z =-a d( P = V2(b + a)F 3 (0,-a). (4.12) 



r-2n 

Rmax — ~ 



Turning now to the finite ip— rod, we note that all solutions we have found have F%/F\ ^ 
1 — 2a/ (b + a) here. Thus, from ( 2.27| ), the coordinate ip possesses a conical excess for 
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< z < b, which is 17 



The proper length of the finite ^-iod is given by 

rb 



L = /^,VMW = ^^/^y^^V^l(M). (4-14) 
Of interest is also the expression for the space-time area spanned by the conical singularity, 



b 



Area = 2vr/? J dz y /f 1 (0,z)h(0,z)f Q (0,z) = 2vr/?y ^— J dz y /F 1 (0,z)F 3 (0,z)F (0,z) . (4.15) 

The GB correction Area\ //3 to the parameter A which enters the thermodynamics can also 
be computed from ( 2.37| ), (f4.1[) , is A = Area/ (3. 

4.2.2 The numerical results 

In the absence of analytical methods to construct EGB black rings, a numerical approach 
of this problem seems to be a reasonable task. We have solved the resulting set of four 
coupled nonlinear elliptic partial differential equations numerically, subject to the above 
boundary conditions. Details on the numerical methods used and on a new coordinate 
system better suited for the numerical study of these solutions are presented in Appendix 
B. 

The problem has two dimensionless parameters, which we have chosen to be b/a and 
a' I a. We note that solutions are equivalent under scaling 

(a,b,a') -> (X 2 a,X 2 b,X 2 a'). 

Mass, area, temperature etc. scale according to their dimensions 

(M, A H , T H , . . . ) - (A 2 M, A 3 Ah, A^Th, . . . ) 

With A 2 = 1/a we obtain the dimensionless quantities 

M = M/a,A H = A R /a 3 / 2 , f H = a l ' 2 T n , . . . 

Starting from black rings of Einstein gravity, we have generated branches of EGB black 
rings by increasing the GB coupling constant a from zero, while keeping the parameters a 
and b fixed. Typical profiles of the solutions are presented in Figures 4 and 5. We note that 
the functions F% are smooth outside of the z— axis, showing no sign of a singular behaviour. 



17 The coordinate if) can of course be rescaled such that its periodicity is 2n on the finite tp— rod. Then 
the interpretation of the solutions is somehow different, since they would describe a static ring sitting on 
the rim of a deficit membrane that extends to infinity (in which case 5 > 0). Since for this choice the 
spacetime is not asymptotically flat, we have prefered to consider the case of a conical singularity localized 
in the bulk. 
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Figure 5. The profiles of the functions Fi used in numerics and of the metric functions /j are shown 
as functions of z for several values of p for a black ring solution with a = 1, b = 2, a' = 0.01. 
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The crucial point here is that the divergent behaviour of the functions ft has already 
been subtracted by the background functions ff. We have verified that the Kretschmann 
scalar stays finite everywhere, in particular at p = 0. 

A number of basic features of these EGB black ring solutions are analogous to those of 
the static black rings of Einstein gravity. In particular, all solutions have a conical excess 
8 on the finite ip— rod (for the choice of Aip = 2n at infinity). Moreover, on the horizon 
the circumference of the S 1 is maximum for z = —a and minimum for z = a. 

The isometric embedding of the horizon is shown in Figure 6 for a family of EGB black 
ring solutions with fixed parameters a = 1 and a' = 0.015 and several values of 6, chosen 
in the interval 1.1 < b < 8. 



Isometric embedding 
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Figure 6. The embedding of the EGB black ring horizon for a family of solutions with a = 1, 
a' = 0.015 and several values of the parameter b. Note, that the conical singularity has been moved 
to the outside of the ring, stretching from its surface to infinity. 

We note, that for this embedding the conical singularity was moved from the finite 
-0-rod to the exterior of the ring, extending from the outer ring circumference (at the 
maximal radius i? m ax) to infinity 18 . Therefore the horizon is smooth and round at the inner 
circumference of the ring, while it exhibits a cone-like edge at the outer circumference at 
the maximal radius i? m ax- As the parameter b is increased, the proper distance from the 
center of the ring to the horizon increases. At the same time the conical excess decreases. 
Consequently, the shape of the horizon becomes more and more spherical with increasing 
b. 

Thus, for a given a', one might hope that the conical excess could be completely 
removed, by going to large enough values of b. This hope, however, is dashed, when the 
parameter space is fully explored. 

Whereas for Einstein black rings the ratio c = a/b exploits the full range < c < 1, 

18 Note that for the choice = 2-k at infinity the embedding would be pseudo-Euclidean. 
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EGB black rings are restricted to the a'-dependent range 

< c mi „(a ) < c < c max (a') < 1. (4.16) 

For a given a' = a' /a, these minimal and maximal values of c can be read from Figure 7. 
This figure, which is one of the central results of this paper, exhibits the maximal value of 
the scaled coupling a!, up to which EGB black rings can be obtained for a fixed value of 
the parameter ratio c = a/b. Therefore it delimits the domain of existence of EGB black 
rings. The coupling a' approaches its maximum value approximately in the middle of the 
interval, i.e. for a/b ~ 1/2. 

As a' — ► a' max , the numerical process fails to converge, although no singular behaviour 
is found there. This result is not a surprise given the black ring - black string connection. 
Heuristically, these EGB black ring solutions may be thought of as being obtained by 
taking a piece of the EGB black string and forming a circle. Thus they would inherit the 
r/j — a' constraint ( |3.18| ) from there. The technical reason which causes the solutions to 
cease to exist at 6t' max is discussed in Appendix B.3. It involves an analytic explanation 
of this fact based on a computation performed in a special coordinate system introduced 
there. Similar to the black string case, the argument in the Appendix uses an analysis of 
the field equations at the event horizon. 

0.05 | 1 1 1 1 1 1 1 1 1 1 




1 ' ' ' ' ' ' ' ' ' 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

a/b 

Figure 7. The domain of existence of EGB black rings is delimited by the maximal value of 
a = a 1 1 a versus the parameter ratio c = a/b. 

In Figure 8 a number of physically relevant properties of the EGB black rings are 
exhibited as functions of the scaled coupling a' for a family of values of the ratio b/a. Also 
included are the curves, showing the values of the respective properties on the boundary 
of the domain of existence of the EGB black rings. 

The first quantity shown is the relative conical excess 5. Only for a' = 0, it covers the 
full range — 1 < 5 < 0. Thus it always differs from zero in the domain of existence of EGB 
black rings, as pointed out above. 

The scaled mass M = M/a is constant for the Einstein black rings, since it depends 
only on a. Therefore the EGB mass curves start for all values of b/a from this Einstein 
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Figure 8. Several physically relevant (scaled) quantities are shown versus a' for several fixed values 
of b/a. The asterisks and crosses indicate the limits b/a — > 1, respectively b/a — ► oo 
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value and increase monotonically until they reach the boundary of the domain of existence 
at the respective maximal value of a' /a. 

The maximal value of the mass, reached on the boundary for b/a — > oo, is finite and 
less than 20% above the Einstein value. Recalling the Smarr-like formula ( 4.10) ) for the 



mass, we conclude that this increase of the mass is basically due to the integral I a i. 

The scaled event horizon area Ah = Ah /a 3 / 2 increases slightly with the coupling a' /a. 
Its boundary line starts from the finite maximal value of the Schwarzschild black hole at 
b/a = 1. Extrapolating the solutions on the boundary to the limit b/a — * oo, we observe 
that the values of Ah of the solutions on the boundary tend to zero as \faTfa. 

The scaled temperature Tjj = a^l^Tn on the other hand exhibits a slight decrease 
with a' /a. (However, the solutions are far away from extremality.) The boundary line for 
the temperature starts from the finite minimal value of the Schwarzschild black hole at 
b/a = 1. Extrapolation this boundary line to the limit b/a = oo shows, that the values of 
Th on the boundary tend to infinity as l/^/a'/a. Consequently, the product of Ah and 



Th, entering the Smarr-like formula ( 4.10 ), stays finite in the limit. 

The inner and outer radii of the ring, i? m i n and -Rmax; exhibit a more pronounced 
dependence on a', in particular for larger values of the ratio b/a. The figure exhibits 
besides the scaled minimal radius i2 m i n also their ratio R m in/ Rm&x- The boundary line of 
the ratio starts from the Schwarzschild black hole value -R m in/-Rmax = and ends at the 
maximally reachable value of R m i n / Rm&x = 1- 

Figure 8 also shows the scaled proper length L of the tp-rod and the scaled thermo- 
dynamical parameter A. For these quantities we observe only a small dependence on a', 
increasing with a' only by a few procent for fixed b/a. A similar picture has been found 



for the quantities A\ and Area\, as given by ( |2.35| ), Q2.37 ), respectively. 



Finally, we would like to return to the intriguing observation, that the conical excess 
5 decreases with increasing the GB coefficient a'. While this result has been unexpected 
(for us), it can be heuristically understood as follows. In the presence of curvature-squared 
terms, the modified Einstein equations leads to an effective stress tensor that involves the 
gravitational field 

= -a'H^ = T^. (4.17) 

Therefore, from some point of view, the quantity a'H\{= —Gfy corresponds to a local 
'effective energy density'. However, this effective stress tensor, thought of as a kind of 
matter distribution, in principle may violate the weak energy condition 19 . This is indeed 
the case, since as one can see in Figure 9, near the horizon this quantity takes negative 
values in some region of the p — z plane. The picture is however quite complicated and 
depends in a nontrivial way on the value of the input parameters a/b and a' . For fixed a' 
and a value of the ratio a/b close to one, the region with a negative 'effective energy density' 
is localized around the topology changing point at p = 0, z = b. This region expands when 
the ratio a/b is decreased. For sufficiently large values of b a more complicated picture 



19 This mechanism has been exploited to construct wormhole solutions in d — 5 EGB theory, see the 
discussion in ill 711 and the references there. 
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emerges, with the occurrence of another region where a'H\ < 0, which is localized around 
the horizon 20 . 

We conclude that, for a black ring, the contribution of the Gauss-Bonnet provides a 
repelling force in addition to that supplied by the conical excess. While this, in principle, 
might have given hope to construct balanced static black rings in EGB theory, these could 
unfortunately not be realized, since the solutions cease to exist before the limit 5 = is 
reached. 

4.2.3 The large b limit and the maximal coupling a' s max 

Let us now address the large b limit and try to understand the character of the solutions 
obtained along the boundary of the domain of existence, when the ratio a/b tends to zero 
along with the coupling a'. 

In Figure 10 we show the ratio fj, = IGttGM/GtvAhTjj versus the ratios R m i n / 'fimax 
and a/b for several values of a'. The dots indicate solutions on the boundary of the domain 
of existence. They represent solutions with the maximal and minimal values for the ratio 
b/a for a fixed value of the coupling a'. 

We note that with decreasing a' the maximal value of a/b increases. Here we expect 
to reach the five dimensional Schwarzschild solution as a' tends to zero and b tends to a. 
In this limit the ratio /i assumes the value fj, = 1. 

On the other hand, the minimal value of a/b decreases with decreasing a', and tends 
to zero as a' tends to zero, whereas the ratio [i tends to a finite value in this limit. In the 
following we show that in this limit the (suitably scaled) EGB uniform black string at its 
maximal value of a' is approached. 

For comparison with the black rings, we have performed a systematic study of the 
EGB uniform black strings. Here we have used a coordinate system which is different from 



( 3.14 ), being related to a limit of the metric ansatz introduced in Appendix B. The line 



element of the black strings is parametrized as 

ds 2 = -f(r)dt 2 + m(r) (dr 2 + r 2 dn 2 2 ) + l(r)d0 2 . (4.18) 

(Note that the d = 4 reduction along the (^-direction gives black holes in isotropic coordi- 
nates.) 

As noted above, the mass and tension of the black strings can be expressed in terms 
of the coefficients entering the asymptotic expansion of the gu and gss components of the 



metric. Thus, from — g tt ~ 1 — , g^ wl + ^,we obtain 



with 



M,/A<t> = M, + ir , T = ^(c, - 2c t ) , (4.19) 



M, = ^ 3 -c . (4.20) 



20 We have noticed a similar picture for EGB black strings. For all solutions there we have found Hi < 
for a region near the horizon. However, the picture is different for the EGB black holes, in which case one 
can prove that Hj is always a positive quantity. A deeper understanding of this difference is still missing. 
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2.5 




Figure 9. Left column: The G\ component of the Einstein tensor near the horizon versus p and z 
for a' — 0.015, a = 1 and b = 1.2, 2.0, 4.0 and 8.0 from top to botton. Right column: The curve 
G\ — versus p and z for the same set of parameters. 
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The Hawking temperature and horizon area of a black string are obtained as 



2vrr 



ro 



(4.21) 



ro 



where ro and A</> denote the isotropic horizon radius and the length of the compact coor- 
dinate, respectively, and the function f(r) = /(r)/(l — ro/r) 2 is finite at the horizon. For 
fixed values of tq, uniform black string solutions exist for all a' s /r^ < a' s max /V 2 ~ 1-44. 

Let us now connect the uniform black strings with the black rings in the limit a/b — > 
0. From the results for the black rings (e.g. Figure 7) we observe that the product 
o'b max (a')/a 2 assumes a finite value when a' /a — > and b m3iX /a — > oo, which corresponds 
precisely to a fraction of the scaled maximal string coupling a^/r 2 ,, i.e. a'6 max (a / )/o 2 — > 

a s,max/^ r 0' 

We therefore introduce the scaling parameter A = ro\^8b/a such that the scaled cou- 
pling a' = A 2 a' = Sbr^a' /a 2 tends to a' s max for a/6 — > 0. The scaled parameters 



A o = r — , 
a 



A 2 /,-r 2 — 
a 2 



then tend to infinity as 6/a — > oo for fixed a. However, the ratio a/v2b = 2ro remains 
finite in this limit. 

In order to obtain the uniform black string limit we also need to scale the angle 
variable <j> for the black rings, since the function also diverges in the limit. Introducing 
cf> = V2b(p = ^rQ(j)b/a and integrating then yields = 8nrob/a, which is the equivalent 
of the asymptotic length of the compact dimension. With this expression we thus find for 
the black rings the scaled Hawking temperature Th, the scaled area per asymptotic length 
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A~u/A(j), and the scaled mass per asymptotic length M/A(f>, 



a 1 



T H = T H /A 



A R /Acl) = \ Aft/Acf) = —\ — Aft/a 

it V a 

M/A4> = \ 2 M/A4> = —M/a , 

7T 




In Figure 11 we show the inverse of the dimensionless Hawking temperature Th = Tutq 
and the dimensionless area An = ^h/AoVq of the black rings on the boundary of their 
domain of existence versus the scaled coupling constant a' = a'/ r o- Also shown are the 
corresponding dimensionless T H = T H ro and A^ = A h /AcVq of the uniform black string 
solutions versus the scaled coupling constant a' s = a' s /rQ. We observe that Th and T H 
assume the same limiting values, when the maximal a! resp. a' s is approached. The same 
holds for the dimensionless entropies An and yl H . We note, that for the black rings the 
values for a' = c/ Simax are extrapolated. The figure also demonstrates, that the product 
j4hTh has only a very slight a'-dependence, i.e. for finite a' it differs only slightly from its 
pure Einstein value (which is 4). 




Figure 11. Left: The inverse of the dimensionless Hawking temperature Th = Th^o and the 
dimensionless area An = Ah/A^q of the black rings on the boundary of their domain of existence 
and the dimensionless Hawking temperature T^ = T H ro as well as the dimensionless area = 
A H /A(/>rQ of the uniform black strings are shown versus a' — a' jr\ and a' s — a' s /r^ : respectively. 
Right: The same for the dimensionless mass M — M / A<j>ro of the black rings and the dimensionless 
mass M* = M*/ro of the black strings. 



To compare the scaled mass of the black rings with the mass of the uniform black strings 
we have to keep in mind that the former are derived for an asymptotically flat space-time 
and therefore have no contribution from the tension. Consequently, the relevant quantity 
of the uniform black strings to compare with is M* in Eq. (^20| ). In Fi gure 11 we also show 
the dimensionless masses M = M / 'AcVq and M* = M^/tq as functions of a 1 /r^ and a' s /rQ, 
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respectively. We observe that both masses assume the same values in the limit where the 
couplings a' s and a 1 tends to the maximal value a' s max . 



4.2.4 The phase diagram 




Figure 12. The dimensionless area An, temperature Th, mass Mh, and the product Ah?h of 
the EGB black rings and Schwarzschild black holes versus a' = a'/r^ as well as the respective 
dimensionless quantities A H Tg, M H and A^T^ of the uniform black strings versus a' s — a' s jr\. 

Let us now compare the properties of the three types of static solutions, the EGB 
Schwarzschild black holes, the EGB black rings and the EGB black strings. To obtain a 
phase diagram for these solutions we exhibit the dimensionless quantities An, Th, and Mh 
of the black holes and black rings together with the corresponding black string quantities 
A s n , T H , and M H in Figure 12. 

We note, that the EGB Schwarzschild black holes connect to the domain of existence 
of the EGB black rings only in the limit of Einstein gravity. Starting from this boundary 
point they form an infinite branch, otherwise disconnected from the black rings' domain of 
existence, with their scaled mass higher and their scaled area lower than the values possible 
for black rings. 
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In terms of the scaled coupling a', the domain of existence of the EGB black ring 
solutions is bounded on one side by the set of black ring solutions with the minimal (re- 
spectively maximal) values of b/a, while on the other side it is bounded by the set of EGB 
black string solutions. The two boundary lines merge and end in a cusp at the maximal 
scaled coupling a' . 

The presence of cusps in the sets of solutions typically signals a change of stability. 
The question of stability of the sets of solutions studied here is, however, beyond our 
current understanding. It would certainly be interesting to see whether the analogue of the 
Gregory-Laflamme instability is present for EGB back strings and leads to associated 
sets of nonuniform EGB back strings. These latter would then be expected to be reflected 
in the emergence of new black rings, which are nonuniform along the S . The presence of 
such solutions would give rise to a much more intricate EGB phase diagram (as would be 
presence of black Saturns, di-rings, etc.). 

Moreover, it would be interesting to consider the phase diagram for dimensions higher 
than five. Construction of the uniform black strings would be straightforward, and one 
might observe a dimension dependence of the maximal GB coupling. Most intriguing would, 
however, be to find a possibility to obtain the corresponding sets of higher dimensional black 
rings. 

5. Further remarks. Conclusions 

In four spacetime dimensions, the Einstein-Hilbert action augmented by a cosmological 
constant term is the unique geometrical action (thus depending only on the metric and the 
curvature tensor) leading to field equations which involve at most second order derivatives 
of the metric. However, this is no longer true if the dimensionality of spacetime is greater 
than four, in which case, from a purely geometrical perspective, there is no compelling 
reason to consider only the Einstein-Hilbert action. 

For d = 5, the most general theory of gravity leading to second order field equations for 
the metric is the EGB theory, which contains quadratic powers of the curvature. Although 
the solutions of this model have been been studied for a long time, most of the literature 
considered only spherically symmetric configurations. Very few EGB exact solutions are 
known and none of them is axially symmetric. 

In this paper we have argued that the approach used to construct d = 5 axially 
symmetric configurations in Einstein gravity within the Weyl formalism, can be used also 
for EGB solutions. This opens the possibility to generalize to EGB theory the known 
solutions with a nonspherical topology of the horizon. 

Although it would clearly be preferable to have analytic solutions (if at all possible), we 
have made some progress in this direction by solving numerically the EGB equations. The 
main purpose of this work was to present a systematic analysis of the static black rings in 
EGB theory. The results of our investigation show that, for a given value of the GB coupling 
constant a', the properties of the solutions are similar to their Einstein gravity counterparts. 
For example, all solutions we have found suffer from conical singularities and such solutions 
are presumably unstable. (The apparent unavoidable existence of conical singularities 
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plague all known asymptotically flat static solutions with a nonspherical topology of the 
horizon (including here also multi-black objects)). Interestingly, the absolute value of the 
conical excess decreases with the GB coupling constant a', but the solutions stop to exist 
for some a' max , before approaching a regular configuration. 

The techniques proposed in this paper can easily be extended to other types of d = 5 
static black objects in EGB theory {e.g. multi-black holes with S 3 topology of the horizon, 
black Saturns, di-rings). Another possible direction to approach in the future is the gen- 
eralization of the results in Section 2 to the case of d — 2 commuting orthogonal Killing 
vector fields, with d > 5. However, these solutions will not be globally asymptotically flat. 
One way to construct d > 5 solutions approaching at infinity the Minkowski background 
and possesing a nonspherical topology of the horizon would be to extend the approach in 
the recent work [37] by including a GB term in the action. 

Let us close this work by briefly mentioning the issue of charged static black rings in 
d = 5 EGB-Maxwell theory. One may hope that the inclusion of some matter fields will 
cure the conical singularity plaguing the vacuum solutions. The natural candidate here 
is a gauge field, in the simplest case an electromagnetic field. Soon after the discovery of 
the vacuum solution, an exact solution describing a U(l) electrically charged static black 
ring was found by several different authors [pSfl , @. However, in Einstein-Maxwell theory, 
the presence of an electric charge alone was found insufficient to stabilize a static black 
ring and prevent it from collapsing, since conical singularities were unavoidable also in the 
charged case. 

A priori, it is not obvious that this result holds also for EGB-Maxwell theory. Therefore 
we have considered a generalization of the EGB black ring solutions discussed in Section 
4 by including a Maxwell term in the action (2.1). Since the Harrison- type generation 
techniques used the construct electrically charged solutions in Einstein-Maxwell theory do 
not hold in the presence of a GB term, one is constrained again to approach this problem 
numerically. Restricting ourselves to a purely electric U(l) potential, A = V(p,z)dt, we 
have performed for the charged case a similar computation to that described in Section 4. 
The boundary conditions satisfied by the metric functions are the same as before, while for 
the electric potential we have imposed ^(0, z) = $ (with $ a constant) on the horizon and 
d p V(p, z)\ p= q = for the rest of the z— axis. At infinity, the electric potential vanishes, the 
electric charge being read from the asymptotic expansion of V(p,z). The solutions were 
found starting with EGB configurations and slowly increasing the value for the electric 
potential on the horizon. 

Although we did not explore yet systematically the full set of parameters, for all 
solutions we have found, the absolute value of the conical excess increases with <£. Thus 
we conclude that the presence of an electric charge alone in very unlikely to stabilize the 
EGB black rings. Nonetheless, we expect that, similar to the Einstein-Maxwell theory, 
the conical singularities would be eliminated by submerging an EGB charged static black 
ring into a background gauge field. A drawback of this construction is that, due to the 
backreaction of the background electromagnetic field, the black ring will no longer be 
asymptotically flat. 

However, the solutions in Section 3 may be viewed as an intermediate step towards the 
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construction of a rotating balanced black ring in EGB theory. In principle, the approach 
in this work can straightforwardly be generalized to the case of spinning solutions pO|j . 
The only obstacle we can see at this moment is the tremendous complexity of the EGB 
equations in the presence of rotation. 
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A. The components of G^ v and 



G v n are: 

A 1 



For the metric ansatz fl2.7p, the essential nonvanishing components of the Einstein tensor 
G p p = e~ 2v ( - v(JJ x + U 2 + Us) + v\U[ + U' 2 + + (if + f/f + C/| (A.l) 
+ (VC/i) • (W 2 ) + (VC/i) • (VC/ 3 ) + (VC/ 2 ) • (VC/ 3 ) + Ux + U 2 + U 3 ) ,(A.2) 



G z z = e~ lv y - v'{U{ + U' 2 + U' 3 ) + v{U x + U 2 + U 3 ) + C/f + U' 2 Z + Ug (A.3) 
+ (VUi) ■ (VC/ 2 ) + (VC/i) • (VC/ 3 ) + (VC7 2 ) • (VC/ 3 ) + + ^2 + c^'l , 



G£ = e~ 2v (v'{U x + U 2 + U 3 ) + u(U{ + U' 2 + U'z) - U' X U X - U' 2 U 2 - U^U 3 -U' x -U' 2 - u£j ,(A.4) 
Gj = e~ 2v ^(VU X ) 2 + (VC/3) 2 + (VC/i) • (VC/ 3 ) + V 2 ^ + V 2 C/ X + V 2 C/ 3 ) , (A.5) 
G% = e~ 2v (^(VC/1) 2 + (VC/ 2 ) 2 + (VC/i) • (VC/ 2 ) + V 2 i/ + V 2 C/i + V 2 C/ 2 ^) , (A.6) 



G| = e- 2v (^(VC/ 2 ) 2 + (VC/3) 2 + (VC/ 2 ) • (VC/ 3 ) + V 2 ^ + V 2 C/ 2 + V 2 C/ 3 j , (A.7) 
The essential nonvanishing components of the tensor are: 

HP = Ae~ Au ( - v'(U x U 3 U' 2 + U 2 U 3 U[ + U 2 U X %) + v{U 3 U[U 2 + U 2 U' X U' S + U X U' 2 U^) (A.8) 

- U X U 2 U 3 {U X + U 2 + U 3 ) + 3(i>U x U 2 U 3 - v'U' x U' 2 U' 3 ) - (plU' x U' 2 + ulu' x u' 3 + UlU[U's) 

- C/i(VC/ 2 ) • (VC/ 3 ) - C/ 2 (VC/i) • (VC/ 3 ) - U 3 (VU 2 ) • (VC/i) 
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HI = 4e" 4iy (i/faUaUli + U 2 U 3 U[ + U 2 UiU' z ) - 0{U 3 U[U' 2 + U 2 UiU^ + UiU 2 U 3 ) (A.9) 

- U{U^(U[ + U' 2 + U' z ) + 3(i/U{U^ - vUiU 2 Uz) - {UgUiU 2 + U?U\U 3 + U{ 2 U 2 U 3 ) 

- U'{{VU 2 ) ■ (W 3 ) - U'i (Wi) • (W 3 ) - t# (W 2 ) • (VC/i)^ , 



flj = 4e~ 4iy f - v'^J 3 U[U' 2 + LW^a + - HhU 3 U' 2 + LW 2 ^ + £>i W 3 ) (A.10) 

- 3(z/cW 2 [/ 3 + vU'^U'^) + (U[Ui + ^)(W 2 ) • (VC/ 3 ) 
+ (E^tr 2 + ^)(VC/i) • (W 3 ) + (t/3^3 + £> 3 )(V*7 2 ) ■ (Vf/i) 



flj = 4e- 4 ^2(Vz,) 2 (V^) • (W 3 ) + z/ (-(U{ + ^)(2cW 3 + C/(^) + U 2 U[ + J7 2 ^) 

+ v (-(f/i + J7 3 )(2^^3 + Ws) + ^3 2 f>i + ^3) - (U[U 3 - U'^f (A.ll) 

- 2(^z> + t/ 3 i/ - U 3 U' Z )U[ - 2{U[v + UW - UiU[)U^ + 2U[U! i 

+ (m 3 - v'u' z ){u'{ - c/i) + (/>[>! - v'u'^u'i - ii 3 ) 

- hu' 2 - u'{u 2 - u 3 u 12 - ugh - hug - u'{u 3 - (vc/i) • (w 3 )v 2 ^ , 



H$ = Ae'^ [2(Vv) 2 (VU 2 ) • (Wi) + u' (-{U{ + U' 2 )(2U 2 h + ^C/() + UlU' 2 + J7 2 2 C/() (A.12) 

+ z> (-(J/! + U 2 ){2U' 2 U[ + U 2 h) + ^ 2 f>2 + C/ 2 2 ^i) " (^2^1 " 

- 2{U[v + UW - UiU[)U' 2 - 2(U 2 v + U 2 v' - U 2 U 2 )U[ + 2U' 2 U[ 

+ (oh - v'u{){u% - u 2 ) + (oh - v'u' 2 )(u'( - h) - hh 2 - ugh - hh 2 - u'{ul 

- u 2 u'( - ugh - (vt/i) • (vu 2 )v 2 i?j , 



Ht = Ae-^ \2(Vv) 2 (VU 3 ) ■ (VU 2 ) + v' (-(U' 2 + Ug)(2U 3 U 2 + U' Z U' 2 ) + U 2 U 3 + U 2 U' 2 ) (A.13) 

+ (-(U 2 + U 3 )(2hU' 2 + U 3 U 2 ) + U?U 3 + U'iU 2 ) - (U!,U 2 - U' 2 U 3 ) 2 

- 2{U' 2 is + U 2 v' - U 2 U' 2 )U' Z ~ 2(U 3 v + U 3 v' - U 3 U 3 )U 2 + 2U' Z U 2 

+ (m 2 - u'u' 2 ){u'i - u 3 ) + (ou 3 - v'u'z){u' 2 ' - u 2 ) - u 3 u' 2 - xj'lij 2 

- U 2 U'i - U%Ul - U 3 U' 2 ' - U 2 U' Z ' - (VL/3) • (VU 2 )V 2 v) , 



where we define 

(VC/) • (Vy) = d p Ud p V + d z Ud z V, V 2 U = d 2 U + d 2 U, (A.14) 
The expressions of these tensors in terms of the functions /j is straightforward. 
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In practice we have solved the following combination of the EGB equations 

EP + E* + E$ + E$-2Ei = 0, EP + E*-2E$ + E$ + E t t =Q, (A.15) 
EP + El + E%- IE* + E\ = 0, EP + El -l(g* + E$ + Ef) = 0, (A.16) 

which diagonalizes the Einstein tensor w.r.t. V 2 U\, V 2 ^, V 2 ^ and respectively V 2 u. 



B. Details on the numerics 
B.l A new coordinate system 

Although we could construct 21 EGB black ring solutions by employing the Weyl-type 



coordinates (p, z), the metric ansatz ( 2.12 ) has a number of disadvantages. For example, 
it has proven difficult to extract with enough accuracy the value of the mass parameter M 
from the asymptotic form of /o and also to study solutions with a/b — > or a/b — ► 1. 

To solve numerically the EGB equations, we have found it more convenient to introduce 



the new coordinates r, 9 and reparametrize the metric ( 2.12 ) as 



ds 2 = -/ (r, 6)dt 2 + T ±- (dr 2 + r 2 d9 2 ) + f^# 2 + / 3 (r, 9) d<p 2 , (B.l) 
r <r < oo , < 9 < | , (B.2) 

where (p, z) are related to (r, 9) by 22 

4 4 4,4 

p= ^sin2#, z= \° cos 29 , (B.3) 

2r 2r 

with r 2 , = a. In these coordinates the horizon is located atr = ro,0<(9< : |. The 
semi-finite and finite ^-rods are mapped to ro < r < oo, 9 = ir/2 and tq < r < rj, 
9 = 0, respectively, while the semi-finite (p-rod is on the interval rj, < r < oo, = 0. Here 
= yjb + \/& 2 — a 2 • 

Next we introduce background functions 

/< = /?4 (B.4) 



analogous to Eq. (4.6). Note that in the coordinates (r, 0) the background functions /q , 
simplify to 



The methods in this case were similar to those used in j37|] to construct d = 6, 7 black holes with 
S 2 x S d ~ i topology of the horizon. 



J Note that this coordinate transformation reduces to (2.15) for ro = 
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Figure 13. The domain of integration for the coordinate system (B.l) is shown for a EGB 
Schwarzschild black hole (left) and a static black ring (right). The relation with the Weyl-type 



coordinates in (2.12) is also presented 



and fi is regular and finite except at the intersection of the ip-iod and the ip-iod. For 
completeness we also include the background functions f® and 



( 1 + (? b 



1 + 



2\ 2 




r b / V r V r / 



-i 



/-K-3 + <i + (^) 4 -(s) 2 — [(? 



where 



12a = -r 



1+ - -2cos26» ^ 



r ^ ) ( 1+ ( rr b) 2C ° S26l (rr b ) ) 



1/2 



The black hole limit corresponds to r& = r , in which case one can easily see that the 
EGB Schwarzschild black hole written in isotropic coordinates is recovered. The relation 
between the coordinates in the new metric form ( [B.lD and the Weyl ones in (2.12) is shown 
in Figure 13 for both black holes and black rings solutions. 

The boundary conditions for the functions F$ follow from the expansions Eqs. ( 2.1S| ), 

oo and the normal 
1 along 9 = 0, 



( 2.19 ) and the assumption of asymptotic flatness. Thus — > 1 as r 
derivatives of all functions vanish on the other boundaries, except -F1-F3 
rf, < r < 00. 



B.2 The numerical methods 



With these parametrisation we solve the resulting set of four coupled non-linear elliptic 
partial differential equations numerically, subject to the above boundary conditions. 



- 39 - 



First, one introduces the new radial variable x = l — r^/r which maps the semi infinite 
region [ro,oo) to the closed region [0,1]. This leads to the following substitutions in the 
differential equations 

r P r — ► (l - x)F <x r 2 F rr — ► (1 - x) 2 F tXX - 2(1 - (B.5) 

for any function F{. 

The equations for F\ are then discretized on a non-equidistant grid in x and 9. Typical 
grids used have sizes 90 x 50, covering the integration region < x < 1 and < 9 < tt/2. 

All numerical calculations are performed by using the programs FIDISOL/CADSOL, 
which uses a Newton-Raphson method. A detailed presentation of the this code is presented 
in |3!J. This code requests the system of nonlinear partial differential equations to be 
written in the form P(x, 9, u, u x ,ug, u x q,u xx ,uqq) = 0, (where u denotes the set of unknown 
functions) subject to a set of boundary conditions on a rectangular domain. The user 
must deliver to FIDISOL/CADSOL the equations, the boundary conditions, the Jacobian 
matrices for the equations and the boundary conditions, and some initial guess functions. 
The numerical procedure works as follows: for an approximate solution vP~\ P(vP->) does 
not vanish. The next step is to consider an improved solution = n^ + An, supposing 
that P(u^' + An) = 0. The expansion in the small parameter An gives in the first order 
= P(nW + An) ps P(u^) + §£(nW)An + g(uW)Au x + ... . This equation can 
be used to determine the correction An^ 1 ) = An. Repeating the calculations iteratively 
( n (3) — n (3) _|_ a u ( 2 ) etc), the approximate solutions will converge, provided the initial 
guess solution is close enough to the exact solution. The iteration stops after i steps if the 
Newton residual P(u^') is smaller than a prescribed tolerance. Therefore it is essential to 
have a good first guess, to start the iteration procedure. 

In each iteration step a correction to the initial guess configuration is computed. The 
maximum of the relative defect decreases by a factor of 20 from one iteration step to 
another. However, for large values of a' convergence is slower. In this case we re-iterate 
the solution until the defect is small enough (about 10 -4 ). Note, that this defect concerns 
the discretized equations. The estimates of the relative error of the solution (truncation 
error) are computed separately. They are of the order 0.001. The errors also depend on 
the order of consistency of the method, i.e. on the order of the discretisation of derivatives. 
For the solutions in this paper, this order was six. We have also monitored the quantities 




(B.6) 



(with Xi, 9i a point of the mesh and en.\ a discretized equation), which provide an average 
error estimate. For most of the solutions, we have found tiha < 10~ n (this holds also 
for the constraint equations E® and E r r — Eg). Most of the errors come from the region 
around the point r = r^, 9 = 0, where the distribution of the points in the mesh should be 
carefully chosen. 

In this scheme, the input parameters are the positions of the rods fixed by a and 
6, (resp. ro and r^) and the value a' of the GB coupling parameter. To obtain EGB 
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black rings, one starts with the Einstein gravity solution as initial guess (i.e. a' = and 
Fi = 1) and increases the value of a' slowly. The iterations converge, and, in principle, 
repeating the procedure one obtains in this way solutions for higher a'. For some of the 
configurations, we interpolate the resulting configurations and use them as a starting guess 
on a finer grid. 

The mass M of the solutions can be determined by extracting the coefficient of the 
1/r 2 decay of the metric component gtt 

-g tt - 1 - c /r 2 , (B.7) 

with Co = 16irGM/3V 3 . Alternatively, the Smarr-like relation ( f4.10| ) between mass, surface 
area, Hawking temperature, and the integral I a i can be employed to determine the mass 
once the other quantities are computed. Comparing this Smarr value for the mass with 
the mass evaluated from the asymptotic decay we find excellent agreement, i. e. deviations 
occur only after 6 digits. 

B.3 The issue of the maximal value of a' 



A = 1-(4 q/p' 



2,1/2 



:1-(4q/p 2 ) 1 ' 2 
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b 1 ? 




b • 3 
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Figure 14. The discriminant A is shown as a function a' for several values of the ratio b/a. 



When a' approaches its maximal value, no singularity shows up. In order to get more 
insight into this behaviour we follow the approach of ref. (2^] and study the expansion of 
the solution near (r = ro, 9 = 0). 

Defining r] = r — r$ we parametrise the solution as 

r 9 2 ~ 

fo = V 2 (h + h , r ri + ho, rr il 2 ) + ( h o,ee + h o , r 0eri + h 0irr o e r] 2 )— , 



/i = (hi + hx,rV + h ltrr r] 2 ) + (hxfie + h 1)r g e rj + h lirree rj 2 )— , 
h = (hz + h 2 , r r] + h 2 ,r r r] 2 )— + (h 2 ,ee + h 2)r ee'n + h^reerj 2 )— , 



h = ( h 3 + h 3 , r r] + h 3)rr r] 2 ) + (h 3) g e + h 3 ^ ee r] + h 3 , rree r] 2 
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where already the boundary conditions at 9 = have been taken into account. Here the 
quantities hi tfll/ ,„ denote constants (in abuse of notation). 

Substitution in the EGB equations yields for the first order expansion coefficients 

ho : r = — W r o , h 0jr ge = —ho,ee/ro , h\ tT = —2h /r , h\fi$ = —ho,ee-r- > 

riQ 

hi ree = — 2/io,06»Aot^ j h 2 . r = , Jl3 >r = , h 3 ^eo = 

For the second order expansion coefficients hi jTr we find a system of quadratic equations. 
After some algebra we find for the coefficients hi irr an equations of the form 

h\ rr +phx,rr + Q = ■ (B.8) 

The coefficients 7io,m h 2 ^ rr and h 3>rr can be expressed in terms of h\ :Tr . 
The discriminant of Eq. (|B. is given by 

p 2 _ [2 ((8h 2 h 3: ee - h 2> eeh3)ho + 2h 2 h z h Qfie ) (2a'h 3j g e - h 3 hir$) a' + h h 2 h 3 l h' 2 r$] h\ 

2 

A real solution to Eq. ( |B.§| ) exists only if ^ — q > 0. We monitored the discriminant and 
observed that the solutions cease to exist exactly when the discriminant becomes negative. 
This is demonstrated in Figure 14. 
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